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Continuous- and discrete-time Glauber dynamics. First- and 
second-order phase transitions in mean-field Potts models 
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Abstract - As is known, at the Gibbs-Boltzmann equihbrium, the mean-field g-state Potts model 
with a ferromagnetic coupling has only a first order phase transition when g > 3, while there is 
no phase transition for an antiferromagnetic coupling. The same equilibrium is asymptotically 
reached when one considers the continuous time evolution according to a Glauber dynamics. 
In this paper we show that, when we consider instead the Potts model evolving according to 
a discrete-time dynamics, the Gibbs-Boltzmann equilibrium is reached only when the coupling 
is ferromagnetic while, when the coupling is anti-ferromagnetic, a period-2 orbit equilibrium is 
reached and a stable second-order phase transition in the Ising mean-field universality class sets 
in for each component of the orbit. We discuss the implications of this scenario in real-world 
problems. 
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Introduction. — The q-state Potts model is one of 
the most important models in statistical physics [I]. For 
example, in the limit g — > 1, the model is equivalent to 
bond-percolation phenomena |5]. More in general, the 
model finds countless applications in computer science, 
like coloring a random graph [2111], and extracting com- 
munities from real-world networks [5] (for a review see 
also [B]). It is perhaps the first candidate model when one 
wants to consider first-order phase transition phenomena 
which, as opposed to second-order transitions, are char- 
acterized by a finite jump of the order parameter. In the 
mean-field version of the model, for example, where any 
spin interact with any other spin, it is easy to show that, 
at the Gibbs-Boltzmann equilibrium, when the number of 
states g > 2, and the coupling is ferromagnetic, a first- 
order transition takes place while, for q = 2, the model 
is equivalent to the Ising model and only a second-order 
phase transition takes place. When the coupling is in- 
stead anti-ferromagnetic, no phase transition sets in, for 
any q, and the system remains in its trivial symmetric solu- 
tion. The same scenario is found when one uses a dynam- 
ical approach like the continuous-time Glauber dynamics 



[5|: after an initial transient time the systems reaches the 
Gibbs-Boltzmann equilibrium, in both the ferromagnetic 
and anti-ferromagnetic case. 

In this paper, after reviewing the mean-field Potts 
model at equilibrium, and along the continuous- 
time Glauber dynamics, we introduce a discrete-time 
"Glauber" dynamics. We then show that, when the 
coupling is ferromagnetic, the differences between the 
continuous- and the discrete-time dynamics are irrelevant. 
However, when the coupling is anti-ferromagnetic, along 
the discrete-time dynamics the system reaches a dynam- 
ical stability by oscillating regularly between two values. 
In other words, when the dynamics is discrete and the 
coupling anti-ferromagnetic, the system never reaches a 
point-like equilibrium, but rather a pcriod-2 stable trajec- 
tory. Moreover, there exists a critical temperature where 
the system undergoes a second-order phase transition in 
the mean- field Ising universality class, for any q > 2. 

Perhaps, the reader familiar with the Ising model {q = 
2) will not find surprising that in the anti-ferromagnetic 
case a discrete-time dynamics can produce stable finite 
oscillations because, at any two dynamic iterations, the 
effective coupling becomes ferromagnetic. Ising models 
with a discrete-time dynamics were studied long ago in the 
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context of 2-state neural networks and cellular automata 
and, later, also in the random field Ising model [T3|. How- 
ever, in the Potts case (g > 2), we find out also the non 
intuitive feature that the nature of the phase transition 
changes dramatically: the fact that the effective coupling 
is ferromagnetic does not return a first-order phase tran- 
sition, but a second-order phase transition. 

Continuous phase transitions in the Potts model with 
q > 3 were observed in complex networks with power law 
exponent 7 < 3 [6], and also on the complete graph in 
the microcanonical ensemble (as opposed to the canoni- 
cal, or Gibbs-Boltzmann, ensemble) [H]- However, the dy- 
namical nature of the continuous transition we present in 
this paper is totally different from the above equilibrium 
cases. We stress also that the discrete-time dynamics we 
shall focus on, is not meant as an approximation of the 
continuous-time dynamics (as is instead usually done for 
practical simulations). There are in fact infinitely many 
remarkable examples where the dynamics is intrinsically 
discrete. Whereas only a continuous-time dynamics can 
represent some description of a system of physical particles 
each other interacting via a physical medium, a discrete- 
time dynamics can represent a system of agents which 
interact via, e.g., exchange of information taking place 
at discrete random times, as in fact occurs in the actual 
world, especially in social or economical contexts, but also 
in ecosystems. As we will show, however, unlike the con- 
tinuous case, in order to be well defined at all times, the 
discrete-time master equation requires imposing a bound 
on how fast free spins can change status (an external pa- 
rameter independent from the Hamiltonian parameters). 

The paper is organized as follows. In the next Section 
we review the mean-field Potts model and discuss the ex- 
istence of an unstable second-order phase transition for 
which we will find later a stable counterpart. In Sec. Ill 
we introduce the continuous- and discrete-time Glauber 
dynamics and discuss in detail their common points and 
differences. We then calculate the critical point and the 
critical behavior of the discrete time-dynamics and com- 
pare with numerical analysis. Conclusions are then draw. 

The mean-field Potts model. We now briefly re- 
view the traditional g-state mean-field Potts model em- 
phasizing some points not often stressed in literature. The 
model is defined through the following Hamiltonian built 
on the fully connected (or complete) graph 



J 



^=-T^Xl'^('^*''^j)' 



(1) 



t<j 



where 5{a, a') is the Kronecker delta function. Let us 
rewrite H as (up to terms negligible for N — >■ oo) 



H 



-E 



1 2 



^<5(cr»,cr) 



(2) 



From Eq. ^ we see that, if J > 0, by introducing q inde- 
pendent Gaussian variables, x^,, cr ~ 1, . . . , g, we evaluate 



the partition function, Z ~ X^icra exp(— /3iJ({(Ti})), as 



(3) 



From Eq. (JS]), for N — >■ cx), by saddle point, we get imme- 
diately the following system of equations 



^PJXc 






,fiJx„ 



a = l,...,q, 



while the free energy density / is given by 



/3/ = -logK^, 



,fiJx, 



E 



l3Jxl 



(4) 



(5) 



to be evaluated in correspondence of the solution of the 
system (j4|). We remind that the term represented by ([5]) 
alone, where the Xo-'s are meant as free variables, is called 
Landau free energy. For each ct, x^ coincides with the 
thermal average J2a 6xp(— /Sff) ^- 5ai,c,/N , i.e., the prob- 
ability to find any spin in the state a. Eqs. ((4]) are sym- 
metric under permutation of the components (xi, . . . ,a;q). 
All the possible solutions of Eqs. (|4]) can be found by 
setting q — 1 components equal to each other and solving 
one single equation. If (zi, . . . , iq) is any permutation of 
(1, . . . , q), then we set 

Xi^ = X, x^. ^y, j = 2,...,q (6) 

where y ~ (I ~ x)/(q — 1), and x satisfies the equation 

1 



1 + (g - l)exp 



9-1 



(7) 



Eqs. ®-(I7]) give rise to a well known phase transition 
scenario [T] : a second-order mean- field Ising phase transi- 
tion sets up only for q = 2, while for any <? > 3 there is a 
first-order phase transition at the critical value (Fig. [T|): 



/3f°-)j: 



2(9-1) 



q 



log{q - 1). 



(8) 



It is important to remind, however, that besides the solu- 
tion giving rise to the first-order phase transition, and the 



trivial symmetric solution xi 



l/q, which, in 



general, can be leading, metastable or unstable, for q > 2 
Eqs. ©-(Ill) have also a non trivial unstable solution giv- 
ing rise to a second order phase transition which sets in 
at the following critical point 



A 



(S.O. unstable) j 



q- 



(9) 



Eq. ([9]) establishes the temperature below which the sym- 
metric solution is no longer metastable or stable. As we 
will see later, Eq. ([9]) has a stable counterpart when J < 
and a discrete-time dynamics is considered. 

In Figs. ([!]) and ([2]) we plot all the possible solutions 
of Eqs. ([7]) and study their nature according to the free 
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energy density ([5]). Note that a solution of Eqs. (0]) is 
leading if it is the absolute minimum for /; metastabl 
if it is a local minimum for /; or unstable when is not 
minimum for / (i.e., at least one eigenvalue of the Hessia: 
of the Landau free energy ^ is not positive) . 



3-state Potts Mean Field 



3-state Potts Mean Field 



ffY?7?Y7?vWw7o.^ — ~ 



Leading 
X,, X Leading 
X| Metastable A 
x^, x^ Metastable A 
X , X, Unstable A 
x^ Unstable A 
X , x^, X Unstable S 
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Fig. 1: (Color online) Magnetizations at equilibrium for the 
case q = 3 and J — 1. The first-order phase transition 
takes place at the critical temperature given by Eq. ((HJ which 
gives Tc ' = 0.36067, while the critical temperature of the 
unstable second-order phase transition is given by Eq. ^, 
which gives Tc ' ' ""'^ ^ cj _ ^^g j^ ^j_^^ figure: "Leading" 
stands for the thermodynamic stable state, i.e., having the 
lower free-energy, "Metastable" stands for the stable state 
having the higher free-energy, while "Unstable S" and "Un- 
stable A" stand for the unstable states having the symmetry 
xi = X2 ~ x-i = 1/3 or not, respectively. The Leading states 
lie on the subspaces Xi^ > Xi^ = Xi^ (in this figure only the 
case xi > X2 = x^ is shown), while the Unstable A states lie 
on the subspaces Xi-^ = Xi,^ > Xi.^ (in this figure only the case 
a;i = 2^2 > xz is shown), where 11,12,13, is any permutation of 
the set of indices {1, 2, 3} and Xi^ -\- Xi^ + Xi^ = 1. Note that 
the asymptotic values of the magnetizations toward T = are 
Xi^ = 1, Xi^ = Xi^ = and Xi-^ = Xi^ = 1/2, Xi^ = 0, for the 
Leading and Unstable A states, respectively. 



Equilibrium versus dynamics and the antiferro- 
magnetic case. — Above wc have analyzed only the 
ferromagnetic case J > 0. For the antiferromagnetic case 
J < 0, it is possible to prove that the system remains 
in the trivial symmetric state [?]• By summarizing, at 
the Gibbs-Boltzmann equilibrium, for q > 3 the system 
undergoes a first-order phase transition only in the ferro- 
magnetic case J > 0, while in the antiferromagnetic case 
J < there is no phase transition. 

A similar situation applies when one considers instead 
a dynamical approach, like the continuous-time Glauber 
dynamics 0. The Glauber dynamics is a single spin-flip 
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Fig. 2: (Color online) Enlargement of Fig. ([T]). Notice, on this 
scale, the region of the metastable A (non symmetric) solution 



in a range of temperatures larger than T^ 



(F.O.) 



0.3607. 



dynamics governed by the following master Eq. 

dp{ai,...,aN]t) _ 
adt 



(10) 



-^^w(crj -^ (T.^)p{<Ti,...,<TN;t) + 

X] X! ^('^i ^ Cri)p(ai, . . . , CTj_i, CTj, CTi+1, . . . CTTv; t), 

where p((7i, . . . , crpf] t) is the probability that the system 
is in the configuration Ci, . . . ,ctjv at time t, w{(j — > a ) 
are transition rate probabilities due to interaction with 
other spins, and a/q is the rate per unit time at which, 
due to the interaction with an heat reservoir, a free Potts 
spin makes transitions from either state to the other 
q — 1. By imposing that, for i — > 00, the system reaches 
equilibrium according to the Gibbs-Boltzmann probability 
p(cri, . . . , UN', t) oc exp(— /37J), with H given by Eq. ([T]), 
and by using the strong law of large numbers, it is easy 
to see that the transition rates must satisfy the following 
functional Eq. (also known as detailed-balance) 



JJx. 






w(a 



BJx 



(11) 



where, as in the previous Section, x^ is the probability to 
find any spin in the state a when i — >• 00. We will keep 
using the same symbol to indicate the probability to find 
any spin ct^ in the state a also for finite t: p{<yi = a]t) = 
Xa{t). Taking into account the normalization condition 
X^ct' '^^'^ — !► (T ) = 1, we will consider the following natural 
solution of Eq. ^ [10] 



w{a 



^PJx 



E." 



RJx 



(12) 



Rates ((T2I) have the remarkable feature, typical of mean- 
field models, that they depend only on the final states. 
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This fact allows us to immediately derive the dynamical 
Eqs. governing the evolution of the Xa{t)'s (also known as 
reduced dynamics 0). By plugging Eqs. (fT2)) in (fTO)) we 
get the following (/-coupled (of which g— 1 arc independent) 
differential Eqs. of the order parameters Xg-{t) 



dXa{t) 

adt 



-Xa{t) 



aPJx„{t) 



Y^^^el3Jx^,(t)- 



(13) 



Eqs. (J13p can be more simply derived also from 
the time-dependent Ginzburg-Landau Eqs. dxa/dt oc 
df{{x'„})/dx„, as it has been done to study the ferro- 
magnetic Potts model in the presence of an oscillating ex- 
ternal field [H]. We here, however, wanted to start from a 
microscopic probabilistic approach. Of course, in the sta- 
tionary limit t — > oo, Eqs. (|13l) are equivalent to the self- 
consistent Eqs. ([4]). Furthermore, among the stationary 
solutions of Eqs. (|13p . the dynamically stable ones coin- 
cide with the minima of the Landau density free-energy at 
equilibrium ([5]). And this holds for both the cases J > 
and J < 0. As a consequence, apart from an initial tran- 
sient, the Potts model at Gibbs-Boltzmann equilibrium, 
or along the continuous Glauber dynamic, do not present 
a significant difference. 

Note that, the Ihs of Eqs. ([T3| and the rate a, that de- 
termine how fast the system reaches equilibrium, arc due 
to the interaction of the spins with the medium environ- 
ment, e.g. heat, while the rhs of Eqs. (|13p stems from the 
interactions of neighboring spins at a given /3. 

In the place of the above continuous-time dynamics, we 
introduce now a discrete-time dynamics taking place at 
the discrete times i = 0, 1, 2, . . . , by the following Eqs. 



Xc{t+ 1) - Xa{t) 



a 



-Xa{t) 



r,PJX„ 



E.'e. 



fiJx^ 



(14) 



Eqs. (|T4)) can be derived from a discrete-time master Eq. 
in close analogy to the continuous case, the difference be- 
ing the assumption that we consider only discrete times 
and then finite difference probabilities. Furthermore, we 
need to work with global transition rate probabilities de- 
fined in the full g^-dimensional space of the N spins. More 
precisely, if we introduce the spin vector a = (ci , . . . , ctat), 
and the associated probability vector p{(t), instead of Eq. 
ITUl). we have 



p{a;t+l)-p{(T;t) 



= -^K'^;i)W^(T^<T')(15) 



-Y,P{^';tW{(T' ^a), 



continuously in time, while in Eqs. (|14p the medium acts 
only at discrete times and, between one time and the next 
one, the spins do not move (alternatively, we can can see 
the spin as frozen). Whereas only Eqs. p^ can rep- 
resent some description of a system of physical particles 
each other interacting via a physical medium, Eqs. (|14p 
can represent a system of agents which interact via, e.g., 
exchange of information taking place at discrete random 
times LJ, as in fact occurs in the actual world, especially in 
social or economical contexts. 



When J > 0, the difference between Eqs. ([T3| and Eqs. 
([T4| is not dramatic as, for t — > oo, their solutions both 
tend to the the same a-independcnt stable point-like sta- 
tionary solutions of Eqs. (j?]) (see Fig. ^). However, 
when J < 0, for i — >■ oo, whereas the solution of Eqs. 
(|13p tends simply to the trivial symmetric one, the solu- 
tion of Eqs. (|14p tends, for t — > oo, to a period-2 stable 
trajectory, i.e., the solution oscillates between two values 
(see Figs. ^ and (O). Furthermore, with respect to 
the temperature, this solution undergoes a second-order 
phase transition belonging to the mean-field Ising univer- 
sality class (as shown in Fig. ([6]), which represents the 
case a = 1). By changing a, the critical behavior remains 
unchanged (and the curves are similar to those shown in 
Fig. ([6])), while the critical temperature turns out to be a 
growing function of a. Analytically, the mechanism can 
be understood and the critical temperature evaluated as 
follows. Let us look again for solutions in the form (j6]). 
From Eqs. (fT4|) we have 



x{t + l) = x(t){l ~a) 



a 



l + (9-l)exp(/?ji^ 



-^Qx(t) 



(16) 



Of course, the stationary solution of Eq. ()16p that we get 
by imposing x{t -I- 1) = x(t) = x, coincides with Eq. ([7]), 
so that, in particular, we already know that such station- 
ary solutions do not undergo any phase transition when 
J < 0. It is worth to remind that, in general, even tough 
a stationary Eq. has a solution, this does not necessary 
mean that the system reaches equilibrium. Eq. (J16p for 
J < represents such a case: there exists the trivial sta- 
tionary solution Xo- = 1/q but, nevertheless, the system 
never reaches equilibrium. Let us consider now the next 
time evolution. From Eq. (|16p we have 



where we have introduced the global transition rates 
in terms of the local weights P^ : W{(t —^ or ) — 
rifci ^(^j ~^ ""i) ^s will be made clear later, for the model 
to be consistent, unlike the continuous case, we must im- 
pose the boundary a < 1. The difference between Eqs. 
(fT3|) and Eqs. ([T4| is that in Eqs. (|13p the medium acts 



^ For simplicity, in Eqs. 1141 1 wc have assumed t = 0,1,2,.... 
Note, however, that the modification to include arbitrary random 
times is elementary since the rhs of Eqs. I I14II involves only an im- 
plicit dependence on the time via the Xi{t)'s. As a consequence, the 
probability distribution for the random times does not affect the be- 
havior of the system. It affects only the times at which the changes 
of the Xi{tys may be observed. 
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xit + 2)=x{t){l-af + 



[l + {q- 1) exp {(3 J 



a{l — a) 
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xit){l- 
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""> ' l + (g l)exp(;9ji--W) 



(g-l)exp(/3ji^ 
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Let us analyze the stationary solutions of Eq. (|17p tha 
we get by imposing x{t -I- 2) = x(t) = x, i.e., the period-i 
stable trajectories. By writing x = 1/g — e, and expandin; 
Eq. pT)) at the first order in e, it is not difficult to se( 
that the solution 6 = becomes unstable when, for givei 
q and a, the temperature satisfies the following Eq. 

a(/3J)2 + 2g(l-a)/3J-g2(2-a) =0. (18 

As expected, when J > 0, the only real root of Eq. (fTSl 
coincides with the critical (unstable) temperature of th( 
equilibrium case Eq. ©. However, when J < 0, the onb 
real root of Eq. P^ is 



/3f-' 



>J = q[l-- 
a 



a < 1. 



(19 



Eq. ()19p says that the critical temperature is a growing 
function of a. Notice that a = 1 returns the same value 
given by Eq. (|9]) with J replaced by | J| but, now, at low 
T, the states live on a subspace that oscillates between 
both the space of the states featured by the stable first- 
order transition that we had for J > 0, and the space of 
the states of the unstable second-order transition that we 
had for J > (compare Figs. H]) and ^). 

Finally, by expanding Eq. (|17p further, up to the thirc 
order in e, and by using Eq. (J19p . we obtain 



Are + {q^ 2)BTe^ + Ce^ + 0(e^ 



0, 



(20 



where r 



m 



(S.O.) 



T)/n 



(S.O.) 



and A, B, and C an 



constants that do not depend on t. From ([2D| it follow; 
that the critical behavior is Ising-like for any a: e oc r^' ■ 
(see enlargement of Fig. ([5])). Notice that in the secont 
term of the Ihs of Eq. (PO]) we made explicit a dependenc( 
on the factor q — 2 so that, the Ising case q = 2, where th( 
second term cancels exactly, is included as a special case 
We come now back to the necessary condition a <1. Ii 
fact, unlike the continuous case, in the Ihs of the maste: 
Eq. p^ we have a difference of probabilities. Whereai 
the derivatives dp/ dt can take any real value, the finit( 
difference of probabilities p(t-|-l)—p(t) that appear in Eqs 
(|15|) . must lie in the interval [-1,1]. Or, in other words, we 
have to ensure that, given that p(t) is a probability, p(t+l) 
is also a probability. Note that, for any real a, both Eqs. 
([To]) or ([TS]) give rise to normalized p(cri, . . . , ctat; i)'s at 
any i, but, in general, it may happen that the p(t)'s can 
take values outside the probability range [0,1]. Let us 
rewrite Eq. (|15|) in compact matrix notation: 

(21) 



From Eq. ([?!]) we see that the matrix appearing to the 
right hand side is a stochastic matrix [i.e. all the entries 
are positive and the rows are normalized to 1) when a < \ 
and, as a consequence the Markov chain p[t) is a well 
defined probability. 
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Fig. 3: (Color online) Continuous- and discrete-time Glauber 
dynamics for a 3-state Potts model with ferromagnetic cou- 
pling J = 1 and a = 1. The continuous and discrete case 
are solutions of Eqs. (|13|l and (|14[) . respectively, with initial 
conditions close to the symmetric solution: a;i(0) = 0.3369, 
3:2(0) = 0.3331, 3:3(0) = 0.33. The chosen temperature is 
T/ J = 1/4, which is lower than the first-order critical temper- 
ature of this case, Tf°-^ /J = 0.3607. 



Glauber dynamics, J=-l, a=l 
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Q \(t) discrete dynamic 

□ X (t) discrete dynamics 

^ X (t) discrete dynamics 

— X (t), x,(t), X (t) continuos dynamics 




p{t + l)^p{t)[aW + {l-a)l]. 



Fig. 4: (Color online) Continuous- and discrete-time Glauber 
dynamics for a 3-state Potts model with anti-ferromagnetic 
coupling ,J = —1 and a — 1. The continuous and discrete 
case are the solutions of Eqs. (|13p and (|14p . respectively, with 
initial conditions as in Fig. [3] The chosen temperature is 
T/J = 1/4, which is lower than the critical temperature of this 
case, Tc ' /\J\ = 1/3, where a second-order phase transi- 
tion takes place. Note that the system never reaches point-like 
equilibrium (see Fig. ((S))). 
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were first-order. We have analyzed here the mean-field 
(j-statc Potts case in detail, but we believe this scenario 
being quite typical and robust with respect to the model 
definition. The results of our analysis shed some light to 
real-world models. Many real-world models, like social 
networks, population growths, economic networks, etc., 
arc in fact characterized by random interactions that take 
place not continuously, but only at certain discrete ran- 
dom times. Furthermore, a portion of these interactions 
are often unfriendly, i.e., anti-ferromagnetic. Therefore, 
at least at the mean-field level, and in the absence of 
quenched disorder, that would involve a glassy behavior 
[T2] , such real- world models are expected to reach a stable 
oscillating behavior and undergo a second-order phase 
transition for each component of the order parameter. 



Fig. 5: (Color online) Enlargement of Fig. 3] Notice, on this 
scale, the oscillating behavior of each component. 
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Fig. 6: (Color online) Magnetizations of the period-2 attrac- 
tor of Eqs. (|14[) for the case q = 3, J — —1 and a — 1. 
Once reached the attractor (for large enough discrete times 
t), each component Xi oscillates periodically taking the two 
values Xi(2t) and Xi{2t + 1) (in this figure only the case 
a;i > ^2 = xq for t even, and xi < a;2 = 2:3 for t odd, are 
shown). Inset: enlargement around the transition point where. 
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we plot also the matching functions y{T) 



for T <T, 

l/3-hl.35(rf ■°->-T)i/2andz(T) = l/3-h0.65(rf ■°->-T)i/^ 

with Tc ' = 1/3, that confirm the Ising-like critical behavior. 

Conclusions. — We have shown that anti- 
ferromagnetic models that evolve along a discrete-time 
dynamics, as opposed to a continuous-time dynamics, 
reach period-2 stable trajectories where, after an initial 
transient, the order parameter oscillates between two 
values each of which undergoes a phase transition. Quite 
unexpectedly, the nature of such a phase transition 
is second-order, even though the only possible phase 
transitions of the Gibbs-Boltzmann equilibrium case 
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